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Abstract 

In the framework of phantom quintessence cosmology, we use the Noether Symmetry Approach to ob- 
tain general exact solutions for the cosmological equations. This result is achieved by the quintessential 
(phantom) potential determined by the existence of the symmetry itself. A comparison between the the- 
oretical model and observations is worked out. In particular, we use type la supernovae and large scale 
structure parameters determined from the 2-degree Field Galaxy Redshift Survey (2dFGRS)and from the 
Wide part of the VIMOS-VLT Deep Survey (VVDS). It turns out that the model is compatible with the 
presently available observational data. Moreover we extend the approach to include radiation. We show 
that it is compatible with data derived from recombination and it seems that quintessence do not affect 
nucleosynthesis results. 

PACS numbers: 04.50.+h, 04.80.Cc, 98.80.-k, ll.25.-w, 95.36,+x 
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I. INTRODUCTION 



Recent analysis of the three year WMAP data yj,|2j,|3j] provides no indication of any significant 
deviations from Gaussianity and adiabaticity of the CMBR power spectrum and therefore suggests 
that the Universe is spatially flat to within the limits of observational accuracy. Further, the com- 
bined analysis of the three-year WMAP data with the Supernova Legacy Survey (SNLS), in yj], 
constrains the equation of state w^, corresponding to almost 74% of dark energy present in the 
currently accelerating Universe, to be very close to that of the cosmological constant value. The 
marginalized best fit values of the equation of state parameter gave -1.14 < w t / e < -0.93 at 68% 
confidence level. Thus, it was realized that a viable cosmological model should admit a dynam- 
ical equation of state that might have crossed the phantom value w = -1, in the recent epoch of 
cosmological evolution. Phantom fluid was first investigated in the current cosmological context 
by Caldwell [4], who also suggested the name referring to the fact that phantom (or ghost) must 
possess negative energy which leads to instabilities on both classical and quantum level |5i |6fl. 
Since it violates the energy conditions, it also could put in doubt the pillars of general relativity 
and cosmology such as: the positive mass theorems, the laws of black hole thermodynamics, the 
cosmic censorship, and causality 171 |8fl. On the other hand, phantom becomes a real challenge 
for the theory, if its support from the supernovae la- Type (SNela) data is really so firm. From 
the theoretical point of view, a release of the assumption of an analytic equation of state which 
relates energy density and pressure and does not lead to energy conditions violation (except for 
the dominant one) may also be useful Q9[|. As for the explanation of the SNela data, phantom is 
also useful in killing the doubled positive pressure contribution in several braneworld models Jioll . 

Phantom type of matter was also implicitly suggested in cosmological models with a particle 
production [11], in higher-order theories of g ravity models Jl2ll . Brans-Dicke models, in non- 
minimally coupled scalar field theories Ill3l . ll4n . in "mirage cosmology" of the braneworld scenario 
II 1511 . and in kinematically-driven quintessence (k-essence) models lllq. Il7ll. for example. Such 
phantom models have well-known problems but, nevertheless, have also been widely studied as 
potential dark energy candidates, and actually the interest in phantom fields has grown vastly 



during the last years and various aspects of pha ntom models have been investigated 118 
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One of the most interesting features of phantom models is that they allow for a Big-Rip (BR) 
curvature singularity, which appears as a result of having the infinite values of the scale factor 
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a(t) — > oo at a finite future. However, as it was already mentioned, the evidence for phantom from 
observations is mainly based on the assumption of the barotropic equation of state which tightly 
constraints the energy density and the pressure. It is puzzling [9] that for Friedmann cosmological 
models, which do not admit an equation of state which links the energy density g and the pressure 
p, a sudden future singularity of pressure may appear. This is a singularity of pressure only, with 
finite energy density which has an interesting analogy with singularities which appear in some 



inhomogeneous models of the Universe [3 
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Recently, phantom cosmologies which lead to a quadratic polynomial in canonical Friedmann 
equation have been investigated Q38Q , showing that interesting dualities exist between phantom 



and ordinary matter models which are similar to dualities in superstrin g co smo 



These dualities were generalized to non-flat and scalar field models 1141 
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53J], brane models Il54ll . and are also related to ekpyrotic models Q55 
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Furthermore, some theoretical studies have been devoted to shed light on phantom dark energy 
within the quantum gravity framework, since, despite the lack of such a theory at present, we 
can still make some attempts to probe the nature of dark energy according to some of its basic 
principles p7|1 . 

Finally, phantom cosmology can provides the opportunity to "connect" the phantom driven 
(low energy meV scale) dark energy phase to the (high energy GUT scale) inflationary era. This 
is possible because the energy density increases in phantom cosmology. Concrete models in this 
sense have been recently elaborated with some interesting results Jssj] . 

In this paper, we want to investigate if the existence of phantom fields can be connected to 
Noether symmetries. Such an issue becomes recently extremely important due to the fact that sev- 
eral phenomenological models have been constructed but, some of them, have no self-consistent 
theoretical foundation. The idea to derive the equation of state from symmetries is not new 115 911 
and recently has been applied to dark energy [60] . From a mathematical point of view, the general 



6l| . Besides, due to 



consideration is that symmetries greatly aid in finding exact solutions [|59L 
the Noether theorem, symmetries are always related to conserved quantities which, in any case, 
can be considered as conserved "charges". Specifically, the form of the self-interacting scalar- field 
potential is "selected" by the existence of a symmetry and then the dynamics can be controlled. 
The equation of state, being related to the form of scalar-field potential, is determined as well. 
However, the symmetry criterion is not the only that can be invoked to discriminate physically 
consistent models but it could be considered a very straightforward one since, as we will see be- 
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low, it allows also to achieve exact solutions. 

In Sect. II, we actually show that phantom fields come out by requiring the existence of Noether 
symmetry to the Lagrangian describing a standard single scalar field quintessential cosmological 
model: we show that it allows a phantom dark energy field, and also provides an explicit form for 
the (phantom) self-interaction potential. Sect. Ill studies how this gives rise to exact and general 
solutions. Also extending the approach to include radiation, we show that it is also compatible 
with the post recombination observational data and that quintessence does not influence the results 
of nucleosynthesis. In Sect IV, we work out a comparison between the theoretical solution and ob- 
servational dataset, as the publicly available data on SNela, the parameters of large scale structure 
determined from the 2-degree Field Galaxy Redshift Survey (2dFGRS)and from the Wide part 
of the VIMOS-VLT Deep Survey (VVDS). In Sect.V, we discuss the presented results and draw 
conclusions. 



II. THE NOETHER SYMMETRY APPROACH 



The Noether Symmetry Approach has revealed a useful tool in order to find out exact solutions, 



in particular in cosmology [61 



62, 



63, 



64]. The existence of the Noether symmetry allows to 
reduce the dynamical system that, in most of cases, results integrable. It is interesting to note that 
the self-interacting potentials of the scalar field Q64Q, the couplings [|6l|l or the overall theory [63], 
if related to a symmetry (i.e. a conserved quantity) have a physical meaning. In this sense, the 
Noether Symmetry Approach is also a physical criterion to select reliable models (see Il63ll for a 
discussion). 

In the present case, let us consider a matter-dominated model in homogeneous and isotropic 
cosmology with signature (-, +, +, +) for the metric, with a single scalar field, 0, minimally cou- 



pled to the gravity. It turns out that the point-like Lagrangian action takes the form 



£ = 3aa 2 - a 3 |ey - V(cf>)j + Da" 3 ^ 



(1) 



where a is the scale factor and the constant D is a constant defined in such a way that the matter 
density p m is expressed as p m = D(a /a) 3y , where 1 < y < 2. For the moment, we will limit our 
analysis to y = 1, corresponding to cosmological dust. The value of the constant e discriminates 
between standard and phantom quintessence fields: in the former case, it is e = 1; in the latter, it 
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is 6 = -1. The effective pressure and energy density of the 0-field are given by 

p<p = e^ 2 -V(<f>), (2) 

p = e l -j> 2 + V(0) . (3) 

P<s> 

These two expressions define an effective equation of state Ws = — , which drives the behavior of 
the model. The field equations are 

2- + H 2 + eh 2 - V(<f>) = 0, (4) 
a 2 

- + 3H<j> + eV'((p) = 0, (5) 
a 

3H 2 =p<f,+ p m , (6) 

where prime denotes derivative with respect to (f>, while dot denotes derivative with respect to time. 

The Noether theorem states that, if there exists a vector field X, for which the Lie derivative of 
a given Lagrangian L vanishes i.e. L X -C = 0, the Lagrangian admits a Noether symmetry and thus 



yields a conserved current H61I1 . In the Lagrangian under consideration, the configuration space 
is M = {a, <p} and the corresponding tangent space is TM = {a, <p, a, (p). Hence the infinitesimal 
generator of the Noether symmetry is 

d d d d 
oa o(p oa d(f> 

where a and fi are both functions of a and and 

da da . 

a = —a + —<p (8) 



da d(f) 

P=£a + £d>. (9) 
da d(p 



The Cartan one-form is 



da dip 
The constant of motion Q = i x 6£ is given by 



6 L = —da + — d<p. (10) 



Q = a(a,<f>)—+f3(a,<f>)- 1 . (11) 
da d<p 
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If we demand the existence of a Noether symmetry, L X -C = 0, we get the following equations, 



-da 
a + 2a— = 
da 

d(p da 

3a + lea — = 

d(p 

3V(<f>)a + aV'(<f>)ap = 



(12) 
(13) 

(14) 
(15) 



We have now to look for conditions on the integrability of this set of equations, limiting ourselves 



to the phantom case (1 .e., e = -1), since the standard case has been already investigated in 1165 
It is possible to assume that a and ft are separable (and non-null), i.e. 

a(a, <p) = AMBiicp), j3(a, 0) = A 2 (a)B 2 (<f>). 



66] 



(16) 



This is not true in general but, in such a case, it is straightforward to achieve a solution for the 
system (I12H15D . It is 



2A cos I \ 



a 



-2V5Asin(i^) 



a 



V((f>) = V sin 



<\ [3 " 



2V2 



(17) 

(18) 
(19) 



which selects the Noether symmetry. 



III. SOLUTIONS FROM NEW COORDINATES AND LAGRANGIAN 

Once that X is found, it is then possible to find a change of variables {a, (p) — » {u, v}, such that 
one of them (say u, for example) is cyclic for the Lagrangian X in Eq. (OQ), and the transformed 
Lagrangian produces a reduced dynamical system which is generally solvable. Solving the system 
of equations ixdu = 1 and i x dv = (where ixdu and i x dv are the contractions between the vector 
field X and the differential forms du and dv, respectively), we obtain: 

a = (v + 9A 2 u 2 )? (20) 

= 2 JI/3 arccos @1) 
Vv + 9A 2 u 2 
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Under this transformation, the Lagrangian takes the suitable form 



£ = D + vV + — + \2A 2 u 2 , 
3v 

where u is the cyclic variable. The conserved current gives 

Q = — = 24A 2 u = B, 
ou 



(22) 



(23) 



which can be trivially integrated to obtain u(t) = Bt + C. We use now the energy condition E £ = 
to find v. We obtain the following differential equation 



v(t) - 3v(0 (D + V v(t) - l2A 2 u(tf) = . 
It is a first order equation which for , i.e. can be factorized into the form 

(p-F l )(p-F 2 ) = 0, 



(24) 



(25) 



being p = v(t) and F t = Fj(t, v). We are then left with solving two first-degree equations p = 
Fi(t, v). Writing the solutions to these first-degree equations as G,(?, v) = the general solution to 
Eq. d24j) is given by the product Gi(t, v)G 2 (t, v) = 0. It turns out that 1 : 

v(t) = £XP( "f^° ( ex P ( V^VW + 48A 2 5 2 V - 4DV f . (26) 
The substitution of the functions a = a(u, v) and (p = <p{u, v) into Eqs. (12012 II) yields 



a(t) 



2 2 (exp - V3VV) (exp ^W t + 4Sco 2 V - 4DV ) 



2\ 



9oj l t l + 



\6V 2 



2 \l- ai'ccos 



3ABt 



9io 2 t 2 + 



(exp-V3VoO( ex P V3V^f+48w 2 V -4DV ) 2 
16VQ 2 



(27) 



(28) 



where we have defined co = AB. Setting a(0) = 0, we can construct a relation among the integra- 
tion constants a), D and V : actually it turns out that D = l+4 *y V ° ■ To determine the integration 
constant co, we set the present time to = 1 . This fixes the time-scale according to the (formally 
unknown) age of the Universe. That is to say that we are using the age of the Universe, to, as a 
unit of time. We then set Oq = a(l) = 1, to obtain 

^8 + ^ - cosh V3Vo 



CO 



V 2 



(29) 



In the following we can set C — in u(t) -Bt + C without loosing generality. 
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The two conditions specified above allow one to express all the basic cosmological parameters in 
terms of Vq, the constant that determines the scale of the potential. It is not directly measurable. 
However we can also strongly constrain its range of variability through its relation with the Hubble 
constant: actually because of our choice of time unit, the expansion rate H(t) is dimensionless, so 
that our Hubble constant H Q = H(t ) is clearly of order 1 and not (numerically) the same as the H Q 
that is usually measured in kms^Mpc -1 . Actually, we can consider the relation 

h = 9.9— , (30) 

T 

where, as usual, h = // /100 and r is the age of the Universe in Gy. We see that H fixes only the 
product hj. In particular, following e.g. yj], we can assume that r = 13.73+°}g, thus we get h < 
0.76 for Hq~\. Since, according to our parameterization, H = —2 , it is 



possible to constrain the range of variability for Vq, starting from H . By means of these choices, 
the exact solutions in Eqs. (T27T) and (|28b can be used to construct all the relevant cosmological 
parameters. In particular 

= -\<Ktf + V(0), (31) 
= ~Htf-V(cp), (32) 

40(o 2 + m 

w<p = . 2 , (33) 

-\(f>{f) +V(<f>) 

«V = ^- 04) 

As it is shown in Fig. (OQ), the model allows an accelerated expansion as indicated from the 
observations, and being a phantom field exhibits a super-quintessential equations of state, with 
w<p < -1 (see Fig. ©), a violation of the weak energy condition (see Fig. ©). Finally, to conclude 
this Section, we present the traditional plot logp^ - log a compared with the matter density (see 
Fig. ©). We see that undergoes a transition from a subdominant phase, during the matter 
dominated era, to a dominant phase. The nowadays accelerated expansion of the Universe can be 
associated to such a transition. Interestingly, both the subdominant and the dominant phases are 
characterized by a constant density behavior. 

Fig.© shows the density parameters of {phantom) quintessence and matter. The fixed point 
regime, characterized by quintessence density parameter equal to unity (Q^ = 1), has not yet been 
reached today. This means that we are living in a transition epoch with Q.j, q ~ 0.7. 
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Figure 1: Redshift dependence of the acceleration a{i): the model allows an accelerated phase of expansion, 
as indicated by the observations. 

A. Including radiation 

A more realistic model can be considered by including also radiation beside dust matter and 
scalar field. In this case, the dynamical equations, as far as we know, do not have analytical 
solutions and it is not possible to analytically reconstruct the Noether symmetry. Due to these 
facts, we will relay on numerical solutions. 

Let us introduce the new independent variable u = log(l + z) = - log(a(t)/a ), where oq is the 
present value of the scale factor. The Einstein scalar field equations can be written in the form 



-2 _ Pm+Pr + V 



(35) 



3 + \^ 



H 2 (f>" = --(Pr+pm) + V (/)' + 



dV 
d(p 



(36) 
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Figure 2: Redshift dependence of the equation of state parameter for the some values of Hq. 

where p r ~ a~ 4 is the energy density of radiation. We numerically solved this system of coupled 
equations, specifying the initial conditions 2 . The results of numerical integration are shown in 
FigsISIZlandE] 

The presence of radiation is hardly changing the behavior of the scalar field, its potential, 
the Hubble constant and the w parameter of the dark energy equation of state. As expected, the 
evolution of the £1 parameters is different. At the initial time (fixed for numerical calculations at 
u = 30), radiation dominates the expansion rate of the Universe, with dark energy and matter being 
sub-dominant. At a redshift of about 5000, the energy density of matter and radiation become 
comparable and, during a relatively short period, the Universe becomes matter dominated. At a 
redshift of about 1 dark energy starts to dominate the expansion rate of the Universe. As result 

2 The initial values are assumed so that u = 30 gives the same value for 0(30), 0(30)' and H(30). 
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Figure 3: Redshift dependence of the equation of state and the function 6(p,p + p$) = p ^ , which allows 
to compare the violation of the energy condition with the super-quintessential expansion. 

(see, Fig. 10), it follows that during the epoch of nucleosynthesis (z ~ 10 9 ) the energy density of 
the scalar field is much smaller than the energy density of radiation. In particular, during such an 
epoch, the kinetic terms in the energy density of scalar field vanishes, and the potential terms is 
constant: in this case, the dark-energy term acts as an effective cosmological constant A, and it 
does not influence the process of primordial nucleosynthesis. 

IV. OBSERVATIONAL DATA AND PREDICTIONS 

The above phantom scalar field model of quintessence, provides an accelerated expansion of the 
Universe which could agree, in principle, to the other cosmological behaviors. To test the viability 
of the model, let us compare now its predictions with some available observational dataset. We 
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Figure 4: Plot of log 10 p^, versus log 10 a (thick line). The dashed lines indicate the log-log plot of log 10 p m oc 
log 10 a" 3 . 

concentrate mainly on different kinds of observational data: the publicly available data on SNela, 
the parameters of large scale structure determined starting from the 2-degree Field Galaxy Redshift 
Survey (2dFGRS)and from the Wide part of the VIMOS-VLT Deep Survey (VVDS). 



A. Constraints from SNela observations 



The model can be constrained by SNela dataset presently available. As a starting point, let 
us take into account the sample of 182 SNIa compiled in Ii67ll . which includes the 21 new SNela 
recently discovered by the Hubble Space Telescope (HST), and combines previous SNela dataset, 



namely the Gold Sample compiled in [68], supplemented by the SNLS dataset II69I1 . 

Following a standard procedure, we perform ax 2 analysis comparing the redshift dependence 



12 



-0.5 0.0 0.5 1.0 

z 



Figure 5: Behavior of density parameters of {phantom) quintessence and matter. The fixed point regime, 
characterized by quintessence density parameter equal to unity (Q^ = 1), has not yet been reached today, 
which can be considered a transition epoch. 

of the theoretical values to the observational estimates of the distance modulus, jj, = m-M, which 
takes the form 



Moreover, the luminosity distance for a general fiat and homogeneous cosmological model can be 
expressed as an integral of the Hubble function as 



where E(z) = -^p is related to the Hubble function expressed in terms of z = ao/a(t) - 1. Let 
us note that the luminosity distance also depends on the Hubble distance c/Hq (which does not 
depend on the choice of the unit of time). Such freedom allows us to fit h or the a priori unknown 



m-M = 5 log D L (z) + 25. 



(37) 




(38) 
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Figure 6: The Hubble parameter as a function of u in the Universe filled in with matter, radiation and scalar 
field. 

age of the Universe t using the SNela dataset. We find \hdl x 2 re d = f° r 1^2 data P omts > an d 
the best fit values are V = 14% which corresponds to H = 0.98 + _° f 4 and Q = 0.68!° °^. We also 
get h = 0.72 + 0.04. In Fig ®, we compare the best-fit curve with the observational dataset. 



B. Dimensionless coordinate distance test 



After having explored the Hubble diagram of SNela, that is the plot of the distance modulus 
as a function of the redshift z, we want here to follow a very similar, but more general approach, 
considering as a cosmological observable the dimensionless coordinate distance defined as : 

f- 



y(z) 



(39) 
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Figure 7: Omega parameters as a function of u in the Universe filled in with matter, radiation and scalar 
field. Q. v is marked in green, Q. r in red and £l m in blue. 



The variable y(z) does not 



depend explicitly on h so that any choice of h does not alter the main 



result. Daly & Djorgovki flVOh have compiled a sample comprising data on y(z) for the 157 SNela 
in the Riess et al. Q68[| Gold dataset and 20 radiogalaxies from [|71|], summarized in Tables 1 and 
2 of Ir70ll . In Ir72ll . they have added the latest SNela data released from the SNLS collaboration 
Ii69l1 thus ending up with a sample comprising 248 measurements of y(z) that we use here. As a 
preliminary step, Daly & Djorgovski have fitted the linear Hubble law to a large set of low redshift 



(z < 0.1) SNela thus finding : 



h = 0.664 + 0.008, 



which is consistent with our fitted value h = 0.72 + 0.04, and with the value H = 72 + 
8 km s _1 Mpc" 1 given by the HST Key project [73] based on the local distance ladder and with 
the estimates coming from the time delay in multiply imaged quasars and the Sunyaev- 
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Figure 8: Plot of log 10 p^ versus log 10 a (thick line). The lines indicate the log-log plot of log 10 p m oc 
log 10 a~ 3 and the log-log plot of log 10 (p r ) oc log 10 (a~ 4 ). 



Zel'dovich effect in X-ray emitting clusters Ir75h . It is interesting to point out that the SHOES 
Team recently completed an extensive new program with the Hubble Space Telescope which 
stream-lined the old distance ladder and observed Cepheids in the near-infrared where they are 
less sensitive to dust. The result was to reduce the total uncertainty in the Hubble constant by 
more than a factor of 2, now to just 4.8 percent uncertainty (h = 74.2 + 3.6) Ir76n . 
To determine the best fit parameters, we define the following merit function : 



1 



N 



— Y 



y(zr, v Q )-yi 



n2 



(40) 



We obtain x 2 red = 1-1 for 248 data points, and the best fit value is V = 14*?, which corresponds to 
H = 0.96^° ^ and £1$ = 0.65 + ° ^ In Fig (flOl) , we compare the best fit curve with the observational 
dataset. Daly & Djorgovski [70] developed a numerical method for a direct determination of the 
expansion and acceleration rates, H(z) and q(z), from the data, using the dimensionless coordinate 
distance y(z), without making any assumptions about the nature or evolution of the dark energy. 
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z 



Figure 9: Observational data from the SNela sample fitted against the model. The solid curve is the best fit 
curve, and the best fits values are Vq - 14+j, which corresponds to Ho = 0.98^^ and = 0.68+™. 

They use the equation 

- q(z) = da/a 2 = 1 + (1 + z) (dy/dz)~ l (d 2 y/dz 2 ) , (41) 

valid for k = 0. Eq. (|4TT) depends only upon the Friedman-Robertson-Walker line element and the 
relation (1 +z) = a /a(t). Thus, this expression for q{z) is valid for any homogeneous and isotropic 
Universe in which (1 + z) = a /a(t), and it is therefore quite general and can be compared with 
any model to account for the acceleration of the Universe. This approach has the advantage to be 
model independent, but it introduces larger errors in the estimation of q(z), since the numerical 
derivation is very sensitive to the size and quality of data. An additional problem is posed by the 
sparse and not complete coverage of the z-range of interest. Measurement errors are propagated in 
the standard way leading to estimated uncertainties of the fitted values. In Fig. (fTTT) . we compare 
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Figure 10: Observational Daly & Djorgovski database ( 17211 ') fitted against the model. The solid curve is the 
best fit curve. 

the q(z) obtained by Daly & Djorgovski from their full dataset with our best fit model. 

C. Growth of density perturbations and observational constraints from galaxies redshift surveys 

A relevant consequence of the presence of a dominant form of dark energy in the Universe, 
in addition to its primary effect on the expansion rate, is to modify the gravitational assembly of 
matter from which the observed large-scale structure originated. In linear perturbation theory, it 



is possible to describe the growth of a generic sma 



through a second-order differential equation [77, 



78] 



1 amplitude density fluctuation 6 M = 8p n Jp n 



5 m + 2H(t)S m - AnGp m 5 m = 0. (42) 
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Figure 11: The allowed region for q(z), obtained by Daly & Djorgovski, from the full dataset (shadow area). 
An approximated polynomial corresponding to a z-window Az = 0.6 is shown with the black-thin solid line. 
With the black-thick dashed lines, it is shown the approximated polynomial fitted to the smoothed data at 
±l<x range, and corresponding to a z-window Az, = 0.4. The red-solid line shows the deceleration function, 
q(z) for our model, corresponding to the the best fit values. 

In Eq. (1421) . the dark energy enters through its influence on the expansion rate Hit). We shall 
consider Eq. (1421) only in the matter dominated era, when the contribution of radiation is really 
negligible. In our model, Eq. (|42l) assumes the form 



2 f-2*(-l - 8V 2 + cosh V3Vo) + V3Vq sinh yJWtf 



+ 



3[ -1 + t 2 + St 2 V 2 - t 2 cosh V3Vo + cosh y/3%t 



(43) 



2 + V (3 + 16V ) - 2 cosh y/W ( 



-8, 



o 



= 0. 



(44) 



71/ 



1 + t 2 + St 2 V 2 - t 2 cosh V3Vo + cosh y/3%t 
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Eq. (1431) does not admit exact analytic solutions. However, since with our choice of normaliza- 
tion the whole history of the Universe is confined in the range t e [0,1], and since we choose 
oj < 1, we can expand the trigonometric functions appearing in Eq. (|43l ) in series around t = 0, ob- 
taining an integrable differential equation, which is a Fuchsian differential equation which admits 
hypergeometric solutions. For the growing mode, we get 



6 + oc t 2 ' 3 2 Fi 



1 7 11 9t 2 Vl 



(45) 



3 6 6 24 + 36V + 192V 2 -24cosh V3Vo_ 
We use such an exact solution to study the behavior of the solution for t - 0, and, mainly, to set 
the initial conditions at t = and numerically integrate Eq. (|43l) in the whole range [0,1]. From 
its solutions, we can define a linear growth rate / that measures how rapidly structures are being 
assembled in the Universe as a function of cosmic time, or, equivalently, of the redshift: 

/-^. (46) 
a In a 

where a is the scale factor. 

The growth index f(z) essentially depends on the value of the mass density parameter at the 
given epoch, Q m . For the cosmological-constant model the dependence is / oc £^); 55 . However, this 
is not valid if the observed acceleration originates from a modification of the equations of General 



aid, 



Relativity; for example, in the Dvali-Gabadadze-Porrati (DGP) braneworld theory 
an extra-dimensional modification of gravity gives f(z) oc Q^ 68 (z). In general, a fitting form 
f(z)f(z) oc Q? m {z) has been shown to be an accurate description for a wide range of models (for 
which flj n (z) itself, not only y, depends on the model). Thus, models with the same expansion his- 
tory H(z) but a different gravity theory, will have a different growth rate evolution f(z) and index 
y. A discrepancy between the measured value of the growth rate and that computed independently 
(assuming General Relativity) from the H(z) yielded by SNela would point out modifications of 
gravity. In Figs. [T2l[T3l we show that, for the parameters of our model, the relation f(z) oc £l^ 68 (z) 
works quite well. Some observational techniques have been suggested to measure f(z) at different 
redshifts. Redshift-space distortions, that is, the imprint of large-scale peculiar velocities on ob- 
served galaxy maps, have not yet been considered in this context. Gravity driven coherent motions 
are in fact a direct consequence of the growth of structure. The anisotropy they induce in the ob- 
served galaxy clustering, when redshifts are used as a measure of galaxy distances, can be quan- 
tified by means of the redshift-space two-point correlation function £(r p ,pi), where r p and n are 
respectively the transverse and line-of-sight components of galaxy separations. The anisotropy of 
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Figure 12: The growth index / in our cosmological model (the solid line), compared with its functional 
dependence from OJ m {z), with y = 0.55 (dashed line). 

f (r p , pi) has a characteristic shape at large r p that depends on the parameter fi = In practice, we 
observe a compression that is proportional to the growth rate, weighted by the factor b L , the linear 
bias parameter of the specific class of galaxies being analyzed. The parameter bL measures how 
closely galaxies trace the mass density field, and is quantified by the ratio of the root-mean-square 
fluctuations in the galaxy and mass distributions on linear scales. A value of /30.49 + 0.09 has 
been measured at z = 0.15 using the 2dF Galaxy Redshift Survey (2dFGRS) sample of 220,000 
galaxies with bias 11821 . 18311 . From the observationally determined {3 and b, it is now straightforward 
to get the value of the growth index at z = 0.15 corresponding to the effective depth of the survey. 
Verde & al. (2001) used the bi-spectrum of 2dFGRS galaxies, and Lahav & al. (2002) combined 
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Figure 13: The relative error between the growth index / in our cosmological model fi„(z), with y = 0.55: 
it turns out that such a relation describes quite well the linear growth rate. 

the 2dFGRS data with CMB data, and they obtained 

b Lverde = 1.04 + 0.11, (47) 
b Llahm = 1.19 ±0.09. (48) 

Using these two values for b the value of the growth index / at z = 0.15 is 

/i = 0.51+0.1, (49) 
f 2 = 0.58 + 0.11 . (50) 

More recently, Guzzo et al. reported a measurement offi = 0.70 + 0.26 at a redshift of 0.77, using 



new spectroscopic data from the Wide part of the VIMOS-VLT Deep Survey (VVDS) Il84h . Using 



a new survey of more than 10,000 faint galaxies, they also measured the anisotropy parameter 
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b L = 0.70 + 0.26, which corresponds to a growth rate of structure at that time of / = 0.91 + 0.36. 
This is consistent with our cosmological model, which gives /(0.77) = 0.97 + 0.12and /(0.15) = 
0.68 + 0.2. However it is also consistent with standard ACDM with low matter density and flat 
geometry, although the error bars are still too large to distinguish among alternative origins for the 
accelerated expansion. This could be achieved with a further factor-of-ten increase in the sampled 
volume at similar redshift. 

V. CONCLUSIONS 

We have investigated the possibility that phantom field dynamics could be derived by the 
Noether Symmetry Approach. The method allows to fix the self-interacting potential of the phan- 
tom field and then to solve exactly the field equations, at least in the case of dark energy and matter 
dominated Universe. The main cosmological parameters can be directly derived starting from the 
general solution. We also worked out a comparison between the theoretical predictions and obser- 
vational dataset, as the publicly available data on SNela and radiogalaxies, the parameters of large 
scale structure determined from the 2-degree Field Galaxy Redshift Survey (2dFGRS)and from 
the Wide part of the VIMOS-VLT Deep Survey (VVDS). It turns out that the model is quite well 
compatible with the presently available observational data. 

Furthermore, we extended the approach to the case including radiation. It can be shown that 
radiation is hardly changing the behavior of the scalar field, its potential, the Hubble constant and 
the w parameter of the dark energy equation of state. As expected, the evolution of the Q. param- 
eters is different. At the initial epochs, radiation dominates the expansion rate of the Universe, 
with dark energy and matter being sub-dominant. At a redshift of about 5000, the energy density 
of matter and radiation become almost equivalent and, for a relatively short period, the Universe 
becomes matter dominated. At a redshift of about 1, dark energy starts to dominate the expansion 
rate of the Universe. It turned out that during the epoch of nucleosynthesis (z ~ 10 9 ) the energy 
density of the scalar field is much smaller than the energy density of radiation, and during such 
an epoch the kinetic terms in the scalar-field energy density vanishes, and the potential terms is 
constant. This meas that the dark-energy term acts as an effective cosmological constant A, and it 
does not influence the process of primordial nucleosynthesis. 

As concluding remark, it is interesting to see that the presence of the Noether symmetry could 
constitute a physical criterion to fix the phantom potential. Such an approach revealed extremely 
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useful also for other classes of models (see, e.g. [63, 



64D). 
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